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Abstract 

We theoretically study the behaviour of a liquid bridge formed between a pair of rigid and parallel plates. 
The plates are smooth, they may either be homogeneous or decorated by circular patches of more hy- 
drophilic domains, and they are generally not identical. We calculate the mechanical equilibrium distance 
of the liquid bridge as a function of liquid volume, contact angle and radius of the chemical domain. We 
show that a liquid bridge can be an equilibrium configuration as long as the sum of the contact angles at the 
two walls is larger than 180°. When comparisons are possible, our results agree well with recent analytical 
and molecular dynamics simulation results. We also derive the effective spring constant of the liquid bridge 
as it is perturbed from its equilibrium distance. The spring constant diverges when the sum of the contact 
angles is 180° and is finite otherwise. The value of the spring constant decreases with increasing contact 
angle and volume, and the rate, at which it decreases, depends strongly on the properties of the two plates. 
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I. INTRODUCTION 



Recently there has been a resurgence of interest in capillary phenomena. One area that has 
come to the fore because of advances in patterning techniques is the behaviour of liquids in contact 
with structured surfaces. These wetting systems are extremely interesting because (i) they often 
show rich static and dynamical behaviours, such as non-spherical drop shapes [1], morphological 
wetting transitions J2, 3], and hysteretic phenomena [4,5], and (ii) further understanding in these 
phenomena leads to smart designs of technological applications. The latter is particularly relevant 
given the current rapid developments of micro fluidic devices fl6|] . 

Here we consider the behaviour of a liquid bridge formed between a pair of solid bodies. In- 
vestigations into such a liquid bridge geometry has a long history. Early studies of this problem 
have been performed, e.g. by Poisson |7|], Goldschmidt and Delauney [|8D in the 19th century, who 
considered surfaces of revolution with constant mean curvature. More recent research in the field 
is motivated by its occurance in many different contexts. Examples include but are not limited 

nn n n 

to capillary condensation and evaporation |9J,LL0|], wet adhesion 111 111 , wet granular matter HI 211 . 



atomic force microscopy II13L Il4ll and low gravity 11151 LL6Q experiments. The existence of these 



liquid bridges often change the systems' physical properties because these bridges exert forces 
onto the surfaces, to which they are attached. 

Due to its broad relevance and implications, a number of aspects of the liquid bridge geometry 
have been addressed. The shape and force displacement relation of the liquid bridge have been 
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shapes dlJ, |20, |21|] and the influence of gravity 1115 



HQ]. 



16, 



The stability of the liquid bridge 



1911 have also been investigated. These 



studies were primarily concerned with chemically homogeneous surfaces. 

Recently, new patterning techniques have been developed that allow the construction of chem- 
ical surface domains on micrometer or even nanometer length scales. The corresponding liquid 
structures are small compared to the capillary length and gravity effects can be safely ignored. 
From the theoretical point of view, bridges between striped surface domains have been consid- 



231] . hi addition, 



ered using both the theory of capillary surfaces as well as lattice gas models |22J, 
polymeric nanobridges between two circular surface domains have also been studied by Molecular 



Dynamics simulations Q24£|. 

In this paper, we will extend these theoretical studies and consider liquid bridges in four differ- 
ent types of slab geometries or slit-pores: (i) First, we will briefly discuss the case of two identical 
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and chemically homogeneous surfaces, which are characterized by the same contact angle; (ii) 
Second, we will consider two chemically homogenous surfaces, which differ, however, in their 
contact angles, a geometry that has been previously studied in Ref. 112511 : (iii) Next, one of the two 
surfaces is taken to contain a circular surface domain, the contact angle of which is smaller than 
the contact angle of the surrounding substrate surface; and (iv) Finally, both surfaces are decorated 
by circular surface domains. 

In the first two cases, the two contact lines of the liquid bridge can move freely along the 
surfaces. In the cases (iii) and (iv), on the other hand, the contact lines may be pinned by the 
boundaries of the chemical surface domains. The pinning 'strength' depends on the wettability 
contrast between the chemical domains and the substrate surfaces. We will focus on a large wetta- 
bility contrast, for which the contact lines are essentially pinned to the domain boundaries. In the 
latter limit, case (iv) of the slab geometry or slit-pore as considered here leads to liquid bridges 
that are rather similar to those formed between two circular plates as studied e.g. in Ref. [21], 
where the contact lines are pinned at the edges of the plates. 

In all four cases, we determine the behavior of the liquid bridges close to their mechanical 
equilibrium states. We define mechanical equilibrium by those liquid bridge configurations that 
exert no force on the solid plates. We will analytically calculate the equilibrium distance of the 
plates and then, as they are slightly perturbed from equilibrium, we shall derive the effective spring 
constant of the liquid bridges. This is the first time, to the best of our knowledge, that the effective 
spring constant has been computed. To verify our analytical results, they will be compared with 
numerical calculations using the Surface Evolver package |26|]. 

We note that this is, in general, a very complicated problem because it depends on the details of 
the surface inhomogeneities which, on real surfaces, will be random in position and spans across 
multiple length scales. Thus, we shall limit ourselves to a subset of these possible scenarios. We 
focus on situations where the liquid bridge solutions still possess an axial symmetry. We shall 
assume that the solids are flat, smooth and parallel to each other, but the solid plates can be either 
chemically homogeneous or decorated by circular patches of more hydrophilic domains, compared 
to the remaining susbtrate surface. 

This paper is organised as follows. In the next section, we will introduce our theoretical de- 
scription as well as the methods that we will use to solve the problem. Then we present our main 
results: the equilibrium distance and the spring constant of the liquid bridge as a function of the 
volume, contact angle, and radius of the chemical domain. A conclusion summarises our results 
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and compares them to relevant studies in the literature. 



H. THEORETICAL DESCRIPTION 

Consider the geometry shown in Fig. [Q A liquid (3 is formed between two parallel plates a 
located at x = —l\ and x = h and is surrounded by a gas phase a. Throughout this paper, we 
shall reserve the subscript 1 and 2 for the left and right walls respectively, as shown in Fig. \T\ 
The plates are smooth, but they can either be chemically homogeneous or patterned with circular 
domains of radius R. In the latter case, let us denote the circular domain as 7 and the surrounding 
domain as 5. We shall further assume that the contact angle of the circular domain 6> 7 is smaller 
than the contact angle of the surroundings 9g. The contact angle is defined as the angle that the 
tangent of the liquid gas interface forms with the solid surface. For the geometry shown in Fig. \T\ 
the contact angles 61 and 9 2 are related to the shape of the liquid gas interface y(x) via 

y'(-h) y'(l 2 ) 

where the prime indicates a derivative with respect to x. 

For a chemically homogeneous plate, the contact angle has a unique value given by the well- 
known Young's equation 

cos# = , (2) 

where the E symbols correspond to the interfacial tensions between the two phases denoted on the 
subscripts. The situation is more complex when the plate is chemically patterned. In this case, the 
value of the contact angle can vary between # 7 and 9s when the contact line is pinned at the domain 
boundary [2]. Inside and outside the circular domain, the contact angles are # 7 and 9$ respectively. 
If we restrict ourselves to axisymmetric solutions, the volume of the bridge is given by 

V = 7T f 2 y 2 (x)dx. (3) 

J—h 

The free energy E of the bridge can be divided into two parts, a bulk term that accounts for the 
liquid-gas interfacial energy and two boundary terms which account for the walls. This decompo- 
sition has the form 



E = 27r£ Q , /3 / dx 

h 

i/(Q 



+ J2 2 W ydypM-VM). (4) 
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FIG. 1 : (a) A schematic diagram of the liquid bridge geometry. The coordinate along the symmetry axis 
is denoted by x, the distance from this axis by y, and the local tilt angle by 4>. The liquid bridge contacts 
two parallel plates at x = —l\ and x = I2 with contact radii y{—l\) and 2/(^2), and contact angles Q\ and #2 
respectively. Note that in general < 9±, 9i < 7r and the contour y(x) does not need to have an extremum. 
(b) Enlarged view close to the contact line as the right plate is displaced from x = I2 (dashed line) to I2 + 5fo 
(full line). The interface shape changes from y{x) to y{x) and the contact line shifts by 5y. 
As before, the £ symbols denote the interfacial tensions, and since we consider axisymmetric 
solutions, the variations of the interfacial tensions Yi a p(y) and T, aa (y) must be axisymmetric too. 
As shown in ||2fl, the contact angle equation © is valid for position dependent interfacial tensions 
as well. 

The pressure difference AP = P a — Pp is a Lagrange multiplier that we need because the 
volume of the liquid bridge is kept constant. As has been shown in the literature (e.g. [ !]), the first 
variation of the free energy functional leads to two equations, an equation for the Laplace pressure 

1 y" 



2M 



(5) 



y(l + y' 2 y/ 2 (1 + y' 2 ) 3 / 2 ' 
and the contact angle equation as given in eq[2l In eq[5l M is defined as the mean curvature of the 
liquid bridge. It is useful to note that for the parameterizations used here, the variations of y(x) 
are automatically tangential to the walls. For general geometries, one has to be careful about the 
variations of the interface shape and of the contact line position. 
Furthermore, using Noether's theorem, we find that the quantity 
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My 2 



(6) 



[l + y>2y/2 

is a constant. In fact, as we shall see below, the constant / is proportional to the force exerted by 
the liquid bridge onto the walls. 

To derive the force, suppose that we have a profile y(x) that satisfies the Euler-Lagrange equa- 
tion (eq[5]) for a given drop volume and separation between the plates. If we now displace the plates 



from [ — Zi, Z 2 ] to [—I1 — 8I1J2+ SI2], the pro file of the bridge changes slightly to y(x) = y(x) + h{x) 
and the energy of the system changes by 
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+ 27r£ a) g 



,/2 



2My 



2/' 



l + y /2)l/2 



(1 + ^)1/2 
2/ ^.,2 



h{x) 
Sy(h) 
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Sh. 



(7) 



x=—li 



[1 + y' 2 ) 1/2 

The first term is always zero because y(x) is an extremal solution. The second and third terms can 
be zero in two cases, when the terms in the bracket go to zero or when 5y = 0. The former is 
satisfied for a smooth homogeneous surface since they correspond to eq[Q while the latter is true 
when the contact line is pinned at the boundary between 7 and 5. In other words, for all the cases 
we shall consider here, only the fourth and fifth terms are important. With this in mind, the force 
can be written as 

5E 



s(i 1 + k) = - 2n ^ f 

-2ttS q/3 ( y{l 2 ) sin 2 - My 2 (l 2/ 



(8) 



-2ttS 



Qfi 



-Zi)sin0i - My 2 (-h 



The force F as given by eq [8] can be understood intuitively as the force acting on one wall and 
is equal to the superposition of the normal force component arising from the liquid-gas interface 
(first term) and the Laplace pressure of the liquid pushing on the wall (second term). When the 
force is positive (negative), the liquid bridge acts to increase (decrease) the separation between the 
two plates. 

Let us now sketch how, in general, we can construct the shape of the liquid gas interface. We 
shall formulate the method in a similar way to Carter lll7n . Since the mean curvature M and the 
force parameter / are constants, we can parametrize the surface using the tilt angle (ft = <t>{x) as 
defined in Fig. [Q which is directly related to the shape contour y = y(x) via 



tan 4> = y' and cos < 



V 



/2\-l/2 



(9) 



Using the latter identity in eq[6l we obtain a quadratic equation for y which has the two solutions 

1 



y± 



2M 



(cos0± (cos 2 0-4M/) 1/2 ). 



(10) 



For small /, these two solutions approach the limiting values y + fa cos <p/M and y_ ~ 0, respec- 
tively. Since a bridge should be characterized by y > for all x, we discard the solution and 
thus obtain the relation 
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y = y+ 



2M 



COS + cos 



AMf) 



between the shape contour y and the tilt angle (f>. 



The x-coordinate can be obtained by using the fact that y'(x) 
and integrating ^ with respect to once, we obtain 



1/2 



dy dcf> 
dd> dx 



sin i 



x 



da- 



cos 2 a 



(11) 



tan <p. Solving for ^ 



(12) 



2M 2M Jo (cos 2 a -AMf) 1 / 2 ' 
Thus far, we have written our solutions in terms of two parameters, / and M . The appropriate 
values for / and M come from our boundary conditions. At the walls, the two tilt angles 



MO = |-02 and 0(-Z 1 ) = 1 -| 



(13) 



which further imply 

cos 62 = sin 02 and cos#i = — sin0i. (14) 

These relations are direct consequences of the geometry in Fig. [T] Substituting these relations to 
eq[T2l we find that the distance between the two walls 



L = h + h, 

COS d\ 



1 



rw/2- 



2M 2M Jo 
cos 9 2 1 r*/2-fc 



da- 
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cos 2 a — AMf) 1 / 2 ' 
cos 2 a 



2M 2M Jo "~ (cos 2 a - 4M/) x / 2 ' 
The second boundary condition comes from the conservation of volume 
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(15) 



(16) 
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III. BRIDGE MORPHOLOGIES FOR VANISHING EXTERNAL FORCES 



The liquid bridge geometry encountered in many physical systems as well as its model ex- 
perimental setups can be categorically divided into two arrangements. In the first category, the 
solid bodies are fixed and the liquid bridge arranges itself to minimize its interfacial energy. In 
the second category, at least one of the solid bodies is able to move and adjust its position. The 
theoretical formulation presented in the previous section is applicable for both categories, but from 
here onwards, we will focus on the latter. 

In this section, we shall calculate the equilibrium distance of the solid surfaces when there is 
a liquid bridge connecting the two plates. We define the equilibrium distance as the distance at 
which the force exerted by the liquid bridge is zero. 

When F — f — 0, the solution for the liquid gas interface simplifies considerably. In this case, 
eq[TT]reads y = cos which implies 

My(l 2 ) = sm(6 2 ) and My(-h) = sin(^) (17) 

at the walls. Since a liquid bridge is characterized by y(l 2 ) > and y(—h) > 0, such a bridge 
must have positive mean curvature M > for contact angles < Q\ < n and < 6 2 < it. 
Furthermore, from eq[[5l we obtain the relation 

cos 9i + cos 9 2 = —ML < 0, (18) 

which is equivalent to 

2 cos (^4^) cos i^r) = ~ ML < • (19) 

Since \9i — 6 2 \ < n and cos ( 6>1 ~ 6 ' 2 ) > for all possible values of the two contact angles, it follows 
from eq [19] that 

9 1 + 9 2 > ix . (20) 

Thus, a liquid bridge can only be in mechanical equilibrium for F — if the sum of the two 
contact angles exceeds 180°. For F — 0, the volume of the liquid bridge also has a relatively 
simple form 

V = vr/^sec 2 6 1 - 1/3) + TrZ^sec 2 6 2 - 1/3). (21) 

Let us now consider the solutions for some special cases. We can categorise them into sym- 
metric (for identical walls) and asymmetric cases (for non identical walls). There are two possible 
scenarios for the former and three for the latter. 
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FIG. 2: Normalised plates separation L/(4V/tt) 1 ^ 3 as a function of contact angle for mechanical equilib- 
rium with external force F = 0. Both solid surfaces are homogenous and have the same contact angle. 
A. Symmetric cases 

When the two walls are identical, 

y{-h) = y{l2) = sm9/M, (22) 
l l = l 2 = L/2 = - cos 9/M, (23) 
V= fL 3 (sec 2 #- 1/3). (24) 

The contact angle is well-defined for homogeneous plates and it is appropriate to characterise the 
equilibrium plate distance as a function of the liquid volume V and the contact angle 9. Fig. 
[2] shows how the plates separation (normalised with respect to the liquid volume) depends on 
the contact angle. The normalised plates separation increases monotonically with respect to the 
contact angle. 

When the solid plates are patterned with circular patches, the contact line may be pinned at the 
boundary between domains 7 and 5. In that case, the contact angle 9 is not unique and its value can 
vary between 6 1 and 9s- On the other hand, the size of the circular patch R sets an intrinsic length 
scale to the system. It is therefore more suitable to express the equilibrium distance as a function 
of V and R. Substituting y(—h) = y(h) = R and using the fact that sec 2 9 = 1 + tan 2 9 = 
1 + 4R 2 /L 2 , we obtain 

V = ^L 3 (2/3 + 4R 2 /L 2 ). (25) 

The equilibrium contact angle itself may be expressed either as a function of the plate separation 
L or the volume V 

tan0 = -2#/L, (26) 
9 
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FIG. 3: Normalised liquid bridge volume V/R? as a function of normalised equilibrium distance L/R. 
Both solid surfaces are patterned with circular domains of radius R. 

V = -27T.R 3 cot 3 0(2/3 + tan 2 9). (27) 

To illustrate how the liquid bridge volume depends on the equilibrium distance of the plates, this 
dependancy is plotted in Fig. [3] 

We have also checked our analytical result against Surface Evolver simulation results (data 
points in Figs. [2] and [3]) and it is clear that the analytical theory and numerical results compare 
favorably. 



B. Asymmetric cases 

Our analytical relations can be used to study cases where the two walls are different, as long 
as the geometry remains axisymmetric. These include (I) two homogeneous plates of different 
contact angles, (II) one homogeneous and one chemically patterned plates, and (III) two chem- 
ically patterned plates. They are more realistic scenarios than the symmetric situations we have 



considered in the previous section. The solutions for all these three cases are given by: 

(/) V = TV ( /3( gec 2 di _ 1/3) + ^3(^2 Q 2 _ 1/3); (28) 

\cos 9 2 J 

L = Z 2 (l + cos0i/cos0 2 ), (29) 

(II) V = -ixll cos 3 6^1 + Rl/l 2 2 fl 2 (sec 2 6 1 - 1/3) + tt/ 3 (2/3 + R 2 2 /l 2 2 ), (30) 
L = l 2 (l - cos0i(l + R 2 2 /l 2 2 ) 1/2 )i (31) 

(III) V = Try^i + l\ - R\ {^-ll + \r\ + \r£) + vr/ 2 3 (2/3 + R\Hl), (32) 
L = l 2 + ^R 2 + l 2 -R 2 . (33) 
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Note that we have parameterised both the volume V and the plates separation L in terms of Z 2 , and 
for case (2), we have assumed that the right wall is chemically patterned with a patch of radius R2 
and the left wall is homogeneous. 

Fig. Ufa) shows how the normalised distance depends on the contact angles for case (I) and 
Fig. HJb-c) how the normalised volume depends on the equilibrium plates separation for case 
(II) and (III) respectively. In Fig. HJa-c), numerical results obtained using Surface Evolver are 
also shown for comparison for several representative cases. The qualitative understandings we 
obtain from the symmetric cases still apply in the asymmetric cases. In particular, the equilibrium 
distance and hence the liquid volume increases monotonically with the contact angle. One new 
and interesting aspect of the asymmetric cases is that for an equilibrium configuration to exist, one 
of the plate can be hydrophilic, or in the case of a patterned plate, one of the contact angles that 
the liquid bridge forms can be less than 90°. This is in agreement with a recent study by De Souza 



et. a/ J25n . However, we recall that for a liquid bridge to be stable, the sum of the contact angles 



6 1 + 6 2 > 180° (eqUB- 

IV. EFFECTIVE SPRING CONSTANTS FOR SMALL EXTERNAL FORCES 

When the plates are displaced from their equilibrium distance, there is a restoring force that 
pushes the system back towards its equilibrium state. This is true as long as one of the plates is 



free to move. One physical example is provided 



by 



an atomic force microscopy (AFM) experiment 



where the position of the AFM tip is not fixed lll3L ll4ll. In this case, one should in principle take 
into account that the AFM tip is curved, with a radius of curvature R tip . However, in the limit of 
y(l) <C i?tip, the curvature of the AFM tip can be neglected. The correction is of order y(l) 2 /Rl ip . 

In the previous section, we found the equilibrium distance for F = f = 0. To obtain the spring 
constant of the liquid bridge, we use a dimensionless expansion parameter e = Mf. Subsituting 
this into eqs [8} [EH EE21 and[T6l and keeping the first order correction terms in e, we obtain 

MF = -27r£ a/3 e, (34) 

My(l) = sin#(l-ecsc 2 fl), (35) 

Ml = - cos 9 - e In (esc 6 + cot 9) , (36) 

and — = J2o=e 1 ,e 2 cos ^ — | cos3 ^ — e cos ^- (3V) 
The barred symbols correspond to the perturbed values of the parameters. 
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FIG. 4: Asymmetric plates, (a) Normalised plates separation L/^V/ir) 1 / 3 as a function of contact angle 
9 for two homogenous plates of different contact angles, (b-c) Normalised volume VjR\ as a function 
of normalised plates separation L/R2 for (b) one homogeneous and one chemically patterned plates, with 
domain radius R2, and (c) two chemically patterned plates, with domain radii Ri and i?2- 
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A. Symmetric cases 



Let us first derive the spring constant for homogeneous and identical plates. We can use eqs 
and [37] to obtain the relations 

. , r 27rcos6* 

AM = -e, (38) 

3VM 2 ' V ' 

_ -AMI - e ln(csc 9 + cot 9) 

" m ' (39) 

for the change in curvature AM and displacement A/ up to first order in e 

Substituting the above equations into eq [34] we find that the force-displacement relation is 
given by 

F = —k 2AI, (40) 

with the spring constant 

k= , a (41) 

2 fel^ + ln (csc£ + cottf)) 

The factor of 2 in front of A/ comes from our convention to define the distance between the plates 

as 21 for the symmetric cases. 

The derivation for the patterned plates is slightly more elaborate. First, we need to solve eqs [331 

and [37] simultaneously to obtain the relations 

(1 + cos 2 9)A9 = 2e cot 6, (42) 
AM_R = cos 9 AO - e esc 9, (43) 

for the change in curvature AM and contact angle AO up to first order in e. Then we substitute 
our results for AM and AO into eq[36]to obtain the displacement Al 

AI -AMI + sin 0A0 - e lnfcsc + cot 9) 

A/ = . (44) 

Finally, after some algebra, we find that the spring constant k in the force-displacement relation is 
given by 

(45) 



In (esc 9 + cot i 



l+cos 2 ( 

We plot in Fig. |5£a) how the spring constant of the liquid bridge depends on the equilibrium 
contact angle for both the homogeneous and patterned surfaces. For a given contact angle, the 
spring constant is stiffer for the patterned surfaces than for the homogeneous surfaces. The spring 
constant also diverges at = 90°, as k oc ^_] T / 2 ) an ^ k oc jqz^jw f° r me homogeneous and 
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FIG. 5: Symmetric plates. Normalised spring constant of the liquid bridge k/2n'E a /s is plotted as a function 
of (a) equilibrium contact angle 6 for both the smooth and patterned surfaces scenarios and (b) normalised 
liquid volume V j R 3 for the patterned surfaces case. 

patterned plates respectively. Then it decreases monotonically with respect to the contact angle. 
This is what one would expect intuitively. The larger the contact angle, the weaker the liquid bridge 
is bounded to the surface and hence it is easier to displace them from equilibrium. Furthermore, 
as we shall see below for the asymmetric cases, a generic feature of our calculations is that the 
spring constant diverges whenever #1 + 6*2 = 180°. Naturally, for two symmetric plates this occurs 
at 9 = 90°. 

For patterned walls, the liquid bridge is better parameterised as a function of volume rather than 
contact angle. The spring constant is plotted in Fig. [5jb) as a function of the normalised volume 
V/R 3 . We have used eq[27]to relate the liquid bridge volume and the contact angle. In terms of 
the liquid bridge volume, the divergence happens when V/R 3 approches zero. 

As in the previous section, we also compare our analytical results against Surface Evolver 
simulations. The latter are shown as data points in Fig. [5] and they agree well with the analytical 
formulas. 
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B. Asymmetric cases 



The same machinery can be applied to cases where the two plates are not identical. One needs 
to expand y, I, M, and 9 to first order in e and then solve eqs [341071 self consistently. Below 
we write down the results for the asymmetric cases. As in the previous section, we consider (I) 
two homogeneous plates of different contact angles, (II) one homogeneous and one chemically 
patterned plates, and (III) two chemically patterned plates. 

mi / /o v w (cos6*i + cos6 ) 2 ) 2 (AK . 

(/) (-/V27r£ Q/3 ) = 3 3^- + (46) 

3 cos 9\ — cos 6*1 + 3 cos 9 2 — cos 6 9 2 

In (esc 9i + cot 9i) + In (esc 9 2 + cot 9 2 ), 

(tt\( b/o^v W (cos9 1 + cos9 2 )cos9 1 cos9 2 

(11) (-k 27rh a p) = — — — + . 2 + (47) 

3 cos 9 1 — cos 6 6*1 + 3 cos 9 2 — cos 6 9 2 sin 9 2 

In (esc 9i + cot 9i) + In (esc 9 2 + cot 9 2 ) — 

( 2 2 sin 3 9 2 cos 9 l \ 



\ sin 9 2 3 cos 9i — cos 3 9 X + 3 cos 9 2 — cos 3 9 2 J 

/ 2 cot 6*2 + 2 cos 9i/ sin 6* 2 + sin 2 9\ cos 9\f sin 6* 2 — sin 6* 2 cos 9\ \ 

\ 2 + 2 cos 2 6*2 + 4 cos 6* 2 cos 6*i + 2 sin 2 6>i cos 9 1 cos 9 2 J ' 

(///) (-k/2n^r = - Am A) - ^g-fM) + + + («) 

sm o\ sm (7 2 sm v\ sm f2 

In (esc 9\ + cot 6 1 !) + In (esc 9 2 + cot 6*2), 

where in the last equation f(a, (3) is defined as 



2 C ot a + 2 ^ + sin2 . /3cos/3 - sin a cos /3 + 2 sin 3 /3 sin a ^""^^ 

f./ q\ sma sin a " " 2H a cos / 



2 + 2 cos 2 a + 4 cos a cos 8 + 2 sin 2 /3 cos /3 cos a + 2 sin 3 /3 sin a p' 3 cos " 

-iXq COS £5 

and the two contact angles are related via R\/ R 2 = sin #1/ sin 6*2. 

The spring constants are plotted as functions of the relevant variables in Fig. [6] for the asym- 
metric cases, with the data points corresponding to Surface Evolver calculations for several repre- 
sentative cases. In Fig. [6£a), it is plotted as a function of the contact angles of the two homogenous 
plates. As we have mentioned in the previous section, for asymmetric plates, equilibrium liquid 
bridge configurations do not require both walls to be hydrophobic. The effective spring constant 
can therefore be defined as long as 6*1 + 6*2 > 180°. When 6*i + 9 2 = 180°, the spring constant 
diverges and as before, its value decreases monotonically with increasing 9\ and 6* 2 . It can be 
shown that the divergence goes as k oc ^ 6i+ l 2 _ w y 

The divergence in the spring constant is also observed for case (II), again when 9± + 9 2 = 180° 
or in terms of volume when the liquid bridge volume approaches zero. This is shown in Fig. 
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FIG. 6: Asymmetric plates. Normalised spring constant of the liquid bridge fc/27r£ a| g is plotted as a 
function of (a) equilibrium contact angle 9 for two different homogeneous plates and (b-c) normalised 
liquid volume V/R^ for (b) one homogeneous and one chemically patterned (domain radius R2) plates and 
(c) two chemically patterned plates with different radii of chemical domains Ri and R2. 
[6tb). In addition, this scenario shows one interesting behaviour. For a given liquid bridge volume, 
there is an optimal contact angle (for the homogeneous plate) at which the spring constant is 
maximum. This is due to the fact that the spring constant depends on the contact angles on both 
surfaces. When 0\ is large, the spring constant is relatively weak as expected. When 9\ is small, 
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the equilibrium contact angle at the patterned wall 6 2 is large and as a result the spring constant 
is again relatively weak. The spring constant is relatively stiff when the values of the two contact 
angles are moderate. 

Finally, the spring constant is shown as a function of the liquid bridge volume and the ratio 
of the radii of the chemical patches for case (III) in Fig. [(Jc). Unlike the two asymmetric cases 
we have considered above, we observe no divergence in this scenario. Instead, the spring con- 
stant reaches a plateau as the liquid bridge volume approaches zero. We previously note that 
the divergence is observed whenever 9\ + 9 2 = 180°. Here the two contact angles are related via 
sin 61 j sin 9 2 — R\/ R 2 and as a result, the condition 9\ + 9 2 = 180° is never satisfied for R 1 ^ R 2 . 
The maximum value of spring constant and the rate at which its value decreases as a function of 
volume depends strongly on the ratio of the radii of the chemical domains. When the ratio ap- 
proaches 1, the maximum spring constant is very large and it decays rapidly with volume. This 
is expected since we ought to recover the symmetric scenario when the ratio is 1 . On the other 
hand, when the ratio is large, the maximum spring constant is small and it decays very slowly with 
volume. To be more precise, here we have set R 2 as the typical length scale and the variation is 
slow relative to this length scale. 



V. SUMMARY 



To conclude, we have theoretically studied the behaviour of a liquid bridge formed between a 
pair of rigid and parallel plates at and close to mechanical equilibrium, corresponding to vanishing 
and small external force F. We have computed the equilibrium distance of the liquid bridge and 
derived its effective spring constant as it is perturbed from equilibrium. The analytical results are 
then verified against Surface Evolver simulations. 



Our results for the equilibrium distance is identical to Carter's H17H for the identically homo- 
geneous plates scenario and compare very well with recent molecular dynamics simulations by 
Yaneva et. al. Il24ll for surfaces patterned with circular domains of radius R. In their paper, Yaneva 
et. al. tentatively used the relation tan 9 = —2R/L derived by Swain and Lipowsky [12211 to 
fit their data for the equilibrium contact angle. Swain and Lipowsky derived this relation when 
they considered a liquid bridge formed between two parallel stripes. This is essentially the two 
dimensional version of the problem considered in this paper. Our results here show that the full 
three dimensional treatment of the problem gives an identical relation and hence explain why the 
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molecular dynamics simulations by Yaneva et. al. fit the theory by Swain and Lipowsky very well 
even though they were done for different geometries. We have also extended our calculations to 
include cases where the two plates are not symmetric. These cases are more realistic and relevant 
to real experimental systems. Whereas in the symmetric cases the plates have to be hydrophobic 
for the liquid bridge to be in mechanical equilibrium with F — 0, this condition is relaxed for the 
asymmetric scenarios as long as the sum of the contact angles is larger than 180°, see eq[T9l 

Our model predicts that surface patterning plays an important role in the effective stiffness of 
the liquid bridge. We find that the spring constant diverges whenever the sum of the contact angles 
is 180°, see Figs. [5|a) and[6£a). This can be achieved for all the cases we have considered here, 
except when the two walls are patterned with different radii of chemical domains. In the later case, 
9\ + 9 2 > 180° and the spring constant is always finite. Another generic property we observe is that 
the spring constant gets softer with increasing contact angle and volume. However, the maximum 
value of the spring constant and the rate at which it decreases depends on the properties of the two 
plates. For example, it is possible to devise a configuration where the spring constant is relatively 
soft for a large range of parameters. 

In our calculations, we have assumed that when the plates are patterned, the contact line is 
pinned at the boundary between the chemical domains, 7 and 5. Such contact line pinning can 
also be provided by topographical patterns, e.g. by circular grooves or posts. Furthermore, the 
contact line may also depin as one varies an appropriate control parameter, such as the liquid 
volume. In the case of chemical domains, the contact line will move outward when 9 > 65 and 
inward when 9 < 9 1 . In this scenario, the results we have derived here for the patterned plates 
are only valid for contact angles 9 that satisfy 9$ > 9 > 7 . When the contact line is completely 
in the 7 and 5 domains, one can simply use the results for the homogeneous plates. One may 
also need to take into account the intrinsic contact angle hysteresis of the surfaces |27|]. Here we 
define intrinsic contact angle hysteresis as the variation in the value of the contact angle due to 
roughness or chemical heterogeneities at length scale which is much smaller than R and V 1 ^ 3 . In 
that case, depending on whether the contact line is advancing or receding, one should appropriately 
substitute the contact angle in our formula with the advancing or receding contact angle. 

Another potential advantage of chemical patterning is that it may offer a simple way to sta- 
bilize liquid bridges against tilt perturbations. For two parallel homogeneous plates, a slight tilt 
or misalignment can result in a substantial displacement of the liquid bridge 12011 . Such structural 
instability is generally an unwanted feature, although in some special cases, such as the transporta- 
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tion of water droplets by shorebirds [12811 . it can be of use. While it remains to be shown by explicit 
calculations, we expect that contact line pinning at the boundaries between the chemical domains 
will hinder the displacement of the liquid bridge, at least for moderate tilt angles. 

The results presented in this paper are relevant for a number of different phy sical systems 
where the liquid bridge geometry plays a role. For examples, Vogel and Steen 111 ill have recently 
proposed a switchable capillary based adhesion device, and our calculations suggest a number 
of ways in which surface patterning can be used to control the stiffness of the capillary bridges 
and hence the strength of the adhesion. In wet granular system 111211 . liquid bridges form and 
they strongly influence the stability of the system. Here the curvatures of the beads may become 
important. Thus, one useful extension of the current work will be to investigate such effects on the 
equilibrium distance and spring constant of the liquid bridges, particularly in situations when they 
are not symmetric. It is also of great interests to investigate the dynamical response of the above 
mentioned physical systems near equilibrium. For such studies, in addition to the results presented 
here, one has to take into account the viscosity of the liquid and the friction of the liquid bridge at 
the parallel planes, as well as solve its hydrodynamic equations of motion. 
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